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DEPARTMENT  OF  THE  ARMY  Rr-.r,.... 

HEADQUARTERS.  UNITED  STATES  ARMY  ARMAMENT/£'?l!lfc^WtT  ^ '•tOO 
ROCK  ISLAND.  ILLINOIS  61201 


REPLr  TO 
attention  OF: 


DRSAR-SA  22  Apr  77 

(SUBJECT:  Communications  Studies  Relative  to  the  Copperhead  (CLGP)  System) 


Mr.  Lee  Waters 

Martin  Marietta  Aerospace  (Orlando  Division) 
M.  P.  - 334 
P.  0.  Box  5837 
Orlando,  FL  32805 


Dear  Mr.  Waters: 

References: 

a.  FONECOiN,  19  Apr  77,  between  Mr.  L.  Waters  (MMA)  and  Mr.  G.  Schlenker 

(DRSAR-SAM),  subject:  Communications  Studies  Relative  to  the  Copperhead 

(CLGP)  Systein. 

b.  ARMCOM  Draft  Technical  Report,  Jan  74,  title:  F‘>eld  Artillery 

Communications  Studie.=  With  Application  to  CLGP  (Inclnsure  1). 

c.  Test  Report  No.  FM255,  Headquarters  MASSTER  (Ft.  Hood,  TX), 

Jail  75,  title:  Forward  Observer  Team  Equipped  With  Ground  Laser  Locator 

Designator  (FOTEGLLD). 

d.  Technical  Report  prepared  by  University  of  Oklahoma  Research 

Institute,  CDOG  Project  65-2,  Feb  66,  title:  Communications  Traffic 

Analysis  of  the  Fire  Support  System  of  Command  and  Control  Information 
System. 


e.  Technical  Report,  US  Army  TRADOC  Systenis  Analysis  Activity, 
Jan  75,  with  Addendum  dated  Mar  75,  title:  Cannon  Launched  Guided 
Projectile  Communications  Study. 

Pursuant  to  the  Ref  a conversation,  DRSAR-SA  transmits  herewith  the 
Ref  b document  as  requested.  It  is  noted  that  this  request  has  the 
concurrence  of  DPCPM-CAHS. 

The  P.ef  fc  study  was  performed  before  the  CLGP  CPEA  was  conducted. 
Study  results  were  presented  to  the  CLGP  Study  Advisory  Group  (S/iG) 


DRSAR-SA 

SUBJECT:  Communications  Studies  Relative  to  the  Copperhead  (CLGP)  System 

and  to  an  interested  PA  staff  group.  The  principal  outcome  of  this 
study  was  an  increased  attention  to  the  commuriications  needs  of  CLGP. 

A direct  manifestation  of  this  improved  awareness  is  the  special 
communications  procedures  used  in  the  CLGP  System  Field  Test,  called 
FOTFGI  I D.  fRef  c).  This  test  shewed  that  by  modifying  the  fire  message 
format  and  commo  procedures  the  present  coriiinumcaciotis  neoworK  waS 
not  overloaded.  Several  factors  are  responsible  for  this  conclusion: 

a.  The  message  arrival  rate  in  FOTEGLLD  was  significantly  smaller 
than  (<l/2)  that  studied  in  Ref  b,  which  is  taken  as  a worst  case  under 
surge  conditions. 

b.  Shortened  message  format  produced  a mean  message  turnaround 
or  service  time  of  about  14  seconds. 

c.  A ten  second  slack  interval  was  permitted  to  alert  the  FO  to 
lase,  as  suggested  by  Ref  b. 

For  your  information  note  Ref  d,  an  early  commo  study  focused  at  the 
battalion  level  with  implications  for  automation.  Also  note  that 
TRASANA  performed  stuaies  relative  to  Copperhead  communi cations  (Ref  e), 
subsequent  to  the  ARMCOM  work. 


Sincerely  yours , 


1 Incl  OTTO  F.  HAASE,  JR. 

as  Acting  Director 

Systems  Analysis  Directorate 

CF: 

Project  Manager,  Cannon  Artillery  V.'eapons  Systems,  ATTN:  DRCPM-CAWS/ 

Mr.  Ed  Manley/Mr.  Ernie  Zimpo/i'ir.  Jim  Williaiiis,  Dover,  NJ  07201 
Director,  US  A^my  Materiel  Systems  Analysis  Activity,  AiTN:  Dt\XSV--G5/ 

Mr.  Dave  Darnnardt,  /'.oerdeen  Proving  Ground,  MD  21005 
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FIELD  ARTILLERY  COMMUNICATIONS  STUDIES 
WITH  APPLICATION  TO  CLGP 


Xnc\  1 


by 

George  Schlenker 

U.S.  Army  Armaments  Command 
Rock  Island,  Illinois 
January  197U 
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SOLUTION  OF  A MESSAGE  QUFJING  PROBLEM 
APPLIED  TO  BATTALION-LEVEL  FIELD  ARTILLERY  FIRE  CONTROL 


Introduction 

Field  artillery  forward  observers  (FOs)  are  presumed 
to  communicate  with  a fire  direction  center  (FDC)  at  ar- 
tillery battalion  level  using  FM  radio.  However,  there 
tti’e , in  gtiiitjx'ax,  iiiox'e  FOo  oiiuii  FM  charuiels  or  discxoto 
frequency  bands.  A typical  situation  is  nine  FOs  and 
three  channels.  The  number  of  FOs  using  the  three  rM 
channels  is  variable,  depending  upon  the  availability  cf 
telephone  lines  and  their  operability  status.  In  the 
following  analysis  the  number  of  users  of  the  commo  net 
is  a discrete  variable  N. 

In  using  the  commo  net  each  of  the  ?'0s  is  regaraed 
as  a customer  competing  for  service  by  an  FM  channel. 

The  service  time,  therefore,  reflects  the  time  required 
by  the  P’O  to  transmit  a fire  message  and  receive  ccn’''^r- 
mation  by  the  FDC.  For  several  reasons  this  service 
time  is  a random  variable  whose  probability  distribu- 
tion is  approximately  gamma  with  shape  parameter  2,  i.e., 
7(2).  A fundamental  property  of  this  distribution  is 
that  it  can  be  considered  as  the  distribution  of  the  sum 
cf  two  identically-distributed  exponential  variables. 

Thus  each  channel  of  service  is  regarded  as  two  stages 
of  service  connected  in  tandem  each  stage  of  which  has 
a service  time  which  is  exponentially  distributed.  The 
service  rate  per  stage  of  service  is  designated  as  a. 
Mince  the  number  of  customers  using  the  commo  net  is 
strictly  finite  and  small,  a customer  usi.ng  or  wanting 
wO  'J.-G  list  nust  be  on  which 

enerates  message  requests.  To  ti-cat  this  • situation. 
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each  FO  is  rep;arded  as  a request  generator  with  arrival 
rate  K.  Therefore  the  total  arrival  rate  for  requests  is 

(N  - n)  \ 

where  n is  the  number  of  FOs  in  the  service  system,  i.e., 
either  using  a channel  or  "in  a queue"  waiting  to  use  a 

cixciri^iO  X • 

The  problem  studied  hero  is  a quantitative  descrip- 
tion of  the  above  queuing  process. 

Definition  of  States  o f the  System 

Define  an  array  of  numbers  such  as  x with  the  follow- 
ing format . 


Stage  of  Service 
Queue  length  Channel  S-j 

X X X 

0^  X X 

C^  X X 

Thus  the  array 

0 10 

0 1 
0 0 

I'epresents:  zero  queue,  one  customer  in  the  first  str.re 
of  service  of  the  first  channel  , one  customer  in  the  se- 
cond stfige  of  the  second  channel,  and  no  customers  using 
the  third  channel. 

V/ith  this  lotTut  the  statf  are  enumerated  ana 
Itablcd  in  Figure  1. 
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Figure  1 

States  of  the  Service  System 


Label 

Configuration 

Customers 
Present  in 
the  System 

0 

0 0 0 
0 0 
0 0 

0 

1.1 

0 10 
0 0 
0 0 

1 

1.2 

0 0 1 
0 0 
0 0 

1 

2.1 

0 10 
1 0 
0 0 

2 

2.2 

0 0 1 
1 0 
0 0 

2 

2.3 

0 0 1 
0 1 
0 0 

2 

3.1 

0 10 
1 0 
1 0 

3 

3.2 

0 0 1 
1 0 
1 0 

3 

3.3 

0 0 1 
0 1 
1 0 

3 

3.U 

0 0 1 
0 1 
0 1 

3 

k.l 

110 
1 0 
1 0 

k 

U.2 

0 0 ] 
1 0 
1 0 

k 

V 


1 

s 

1 

i I 
< i 


i 


1 


i 

i 

i 


I 

i 

i 

f 

i 

j 
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Figure  1 
(Continued) 


Label 

Configuration 

Customers 
Present  in 
the  System 

1;.3 

10  1 
0 1 
1 0 

k-k 

10  1 
0 1 
0 1 

k 

In 

general , 

for  j customers  present 

in  the  system, 

3 ^ J 5 

N, 

j-1 

j-3  1 0 

1 0 
1 0 

j 

j-2 

0-3  0 1 

1 0 
1 0 

j 

J.3 

j-3  0 1 

0 1 
1 0 

0 

j-U 

j-3  0 1 
0 1 
0 1 

j 
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Changes  or  transitions  between  these  states  occur 
stochastically.  The  conditional  probabilities  for  the 
occurance  of  these  transitions  in  a differential  time 
interval  dt  can  be  represented  in  a State  Transition 
Diagram  shovm  in  Figure  2.  The  direction  taken  oy  the 
transition  is  indicated  by  the  direction  of  the  arrow 
(arc)  between  labeled  nodes,  each  of  which  represents 

C.  o ^ PT*  T>  ^ P P T*  53  t3^P 

transitions  are  placed  adjacent  to  the  proper  arcs. 
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state  Equations 

The  set  of  aifferential  equations  describing  the  dy- 
namics of  this  service  system  can  be  written  by  inspect- 
ing the  state  transition  diagram.  In  writing  the  dif- 
ferential equation  for  the  probability  of  being  in  state 

n at  time  t,  P (t),  one  examines  arcs  at  the  n th  node, 
n 

One  simplyf  assigns  the  rate  parameter  of  a particular 
arc  entering  tnat  node  as  a positive  coefficient  to  tne 
probability  associated  with  the  origin  of  that  arc.  For 
arcs  leaving  the  n th  node,  a negative  value  is  assigned 
the  rate  paraneter(s)  which  is  (are)  used  as  coeff icient( s ) 
of  P (t).  Thruout  the  netv;ork,  rate  parameters  are 
treated  as  proportionality  constants  of  flow  rates  and 
Kirchoff's  law  is  satisfied.  As  a consequence  of  satis- 
fying Kirchoff's  law,  the  sum  of  all  elements  in  each 
columji  of  the  coefficient  matrix  of  the  differential 
equations  must  equal  zero.  The  last  property  can  be 
used  to  check  the  validity  of  the  aiff erenaj al  equations. 

In  the  follow;' ng  equations  v;e  omit  the  period  in 
the  state  labels  used  in  Figure  1.  Additionally  the 
functional  dependence  on  time  of  the  probabilities  has 
been  suppressed. 


P = - N P + u P-,  „ 
o o ' 12 


= N V P^  - [(N-1)  V + gj  P^^  + a P^p 


Pl2  = A - [(N-1)  ^ ...  .J  . 2 a Pp3 


PpT  = (iJ-1)  P^^  - t(N-2)  \ + 2 ij  Pp^  + a P^p 


Pp2  = (N-1)  P^p  + 2 i-  Pp. 


[ (N-2)  X 2 ij  ppp  2 J.  P33 
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^23  = ^22  ■ ^ ?23  + P3i^ 

= (N-2)  \ P21  - [(N-3)  \ + 3 u]  P^^  + a P^2 
P32  = 3 u-  P31  ^ (N-2)  A P22  - [(N-3)  A + 3 P 

+ 2 p,  P^^ 

P^^  = 2 ti  P^2  + ^^-2)  ^ ^23  " \ + 3 P-]  P 

- 3 - 

P3,  = H ?3J  - [(n-3)  + 3 n) 

For  a nvimber  in  the  system  n with 
14.  $ n 5 N - 1 

^lOn+1  “ (N-n+1)  A P]_o(n-l)  + l 

- t(N-n)  \ + 3 Pion+1 

■*■  ^10(n+l)  + 2 

^lOn+2  = ^ ^lOn+1  (N-n+1)  A Pio(n-l)  + 2 

- [(N-n)  \ + 3 + ] Pion+2 
+ 2 .1  P 


10(n+l)+3 


^lOn+3  = ^ ^ ^lOn+2  (N-n+1)  \ Pio(n-l)+3 

- t(N-n)  W 3 a] 

^ ^10(n+l)+4 

^’lOn+4  ■ ^lOn+3  ^ ^10(n-l)+i|. 

- [(N-n)  X . 3 a]  Pion.4 

The  lapt  four  states  of  the  system  are  described  by 


P - V 

^1011+ 1 - ^ 

^10(K-1)+1  ■ 

^ ^lON+1 

• 

^lOH+2  = ^ 

^lON+1  ^ 

^10(N-l)+2  ' 

3 '-t 

^lON+2 

^1011+3  * ^ 

^’lON+2 

^10(N-l)+3  ~ 

3 ^ 

^lON+3 

^lON+4  " ^ ^lGIJ+3  ■*■  " ^10(N-l)+4  " 3 ^:ON+lj. 
Niimerical  Formalism 

In  using  a computer  program  to  solve  the  above 
differential  equations,  it  is  necessary  to  place  the 
state  variables  in  a sequential  array  U(i,j),  where  1 
is  indexed  from  one  to  the  number  of  states  describinE< 
the  system  and  where  j r I for  the  state  probability  and 
j ■ 2 for  the  first  derivative  of  the  state  probability. 

The  number  of  states  for  the  system,  i.e.,  the 
maximum  vnl\je  of  the  index  i above,  depends  upon  the 
number  of  customers  N as  follows. 
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max(  i ) = U-  N - 2 


(2) 


In  the  computer  program  developed  by  Stuart  Olson  (WECOM 
Report  No.  PAA-TRI-72) , [1]  an  array  is  required  for 
work  space  W(k).  The  dimension  of  k in  this  array  is 
given  as 

ma^rlk)  = 3 max(i)/2 

max(k)  = 3 (2  N - 1)  (3) 

An  association  between  the  mathematical  notation 

Pl3(t) 

and  the  computer-oriented  notation 
U(k,l) 

is  provided  below. 

U(l,l)  - Pq 
U(2,l)  = P;l;l 
U(3,l)  - 
U(U,1)  = P^i 
U(5,l)  = ?22 
U(6,l)  = P23 
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Generally,  for  k 5 7> 


k = i4.i+J-6 


so  that 


Li 


UO;  i + j - 6,  1)  = , 3 € i $ N,  (4) 

1 ^ j « I4.. 

Collections  of  States 

As  implied  in  Figure  1,  certain  states  may  be  treated 
collectively  to  describe  a measurable  property  of  the 
system.  For  example,  the  probability  of  finding  n cus- 
tomers present  in  the  service  system  is  obtained  as 
follows 


P{n  in  the  system}  ■ p 


with 


Pq  = U(l,l) 


= U(2,l)  + U(3,l) 

P2  = U(l;,l)  + U(5.1)  + U(6,l) 


Pn  = 2 U(j+n+3,l) 

,1  = 1 


3 s n ^ N . 


The  number  in  the  queue,  i.e,,  waiting. to  capture 
a channel,  is  just  three  less  than  the  n\imber  in  the 
system  when  this  is  three  or  greater.  The  expected  num- 
ber in  the  queue 

E[QJ  =2^  (n-3)  P„  (6) 

n=l|. 


Parameter  Values 


In  many  tactical  situations  the  message  arrival 
rate  within  the  commo  system  will  be  time-dependent.  It 
is  nevertheless  instructive  to  examine  system  dynamics 
in  response  to  a sudden  initial  arrival  rate,  K,  held 
constant  thereafter.  The  system  is  assumed  to  be  vacant 
initially.  All  results  presented  in  the  following 
graphs  were  ootained  for  a,  = O.Olbo?  sec  . The  corres- 
ponding mean  time  between  arrivals  per  PO  (MTBA)  of  60 
seconds  is  regarded  as  the  shortest  interval  of  interest. 

An  interesting  range  of  service  rates  for  the  CLGF 
system  is  from  0.04  'to  0.10  sec"  . These  correspond  to 
mean  times  to  serve  (MTTS)  of  50  and  20  sec.,  respec- 
tively [2]  . 

Results  and  Interpretation 

Using  expected  queue  length,  P[Q],  as  a measure  of 
the  state  of  the  service  system,  one  can  compare  the 
approach  toward  stationary  queuing  by  inspecting  a plot 
of  E[ Qj  versus  time.  This  type  of  function  is  displayed 
in  Figure  3 for  the  case  of  9 POs  and  3 channels.  Note 
that  equilibrium  queuing  develops  substantially  earlier 
for  a MTTS  of  20  sec  than  for  a MTTS  of  50  sec.  The 
stationary,  expected  queue  length  is  also  shown  to  be 
quite  sensitive  to  the  service  rate. 

Another  description  of  the  system  which  exhibits 
interesting  transients  is  the  probability  that  there  are 
n customers  in  the  system  at  time  t.  This  is  shown  in 
Figure  4 for  the  case  of  5 POs  and  3 channels  with  a 
MTTS  of  50  sec.  Also  shown  in  this  figure  is  the  expected 
queue  length.  By  comparing  Pigui-es  3 and  4.  one  can 
appreciate  the  sensitivity  of  E[Q]  to  the  number  of  customers. 


For  a system  having  5 FOs  and  3 channels,  the  effect  of 
reducing  the  MTTS  to  20  sec  is  shown  in  Figures  5 and  6. 
The  equilibrium  probability  that  four  customers  are  pre- 
sent is  less  than  0.02  whereas  with  a MTTS  of  SO  sec, 
this  probability  is  O.II4.  (from  Figure  4). 

For  a battalion-level  commo  system  with  all  nine 
FOs  using  the  three  FM  channels  with  equal  intensity, 
tne  stationary  prooaoixity  tnat  an  rO  must  wait  for  a 
channel  and  the  expected  queue  length  are  shown  in 
Figure  7 as  functions  of  the  mean  time  to  serve.  Clear- 
ly if  such  a system  were  to  be  feasible  with  the 
Cannon-Launched  Guided  Projectile  System  (CLGF),  the 
MTTS  would  have  to  be  very  short  --  less  than  12 
seconds  to  assure  a 90%  chance  of  capturing  a channel 
promptly.  If  the  channels  were  divided  so  that  three 
FOs  were  compelled  to  share  one  channel,  the  resulting 
queue  behavior  would  differ  from  the  9 FO  - 3 channel 
case  evcii  oi.o  thx  i*ati.o  ox  rOs  to  criannexs  i-s  identical. 
Results  for  the  3 MO  - 1 channel  case  are  shovm  in 
F’igure  8.  These  results  should  be  compared  w_th  those 
of  Figure  7.  Other  thinfts  being  equal,  the  segregation 
of  channels  is  seen  to  increase  the  expected  queue  and 
probability  of  waiting.  As  shown  in  Figure  9,  the  ex- 
pected time  to  wait,  EI.W],  and  the  expected  time  to 
wait,  given  a requirement  to  wait,  are  both  strong  func- 
tions of  the  mean  time  to  serve  for  the  case  of  3 MOs 
and  1 channel.  The  equations  used  in  developing  the 
results  for  the  case  of  3 MOs  and  1 channel  are  shown 
in  the  next  section  (p.  22). 
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Sepire^ated  Communicat ions 

In  this  section  we  consider  the  situation  in  which 
the  conventional  three  Fli  channels  per  field  artillery 
battalion  are  allocated  among  the  FO  users  (customers) 
so  that  each  channel  is  given  to  a separate  set  of 
users  which  pre-empt  all  others  in  using  that  channel. 
Assume  these  customers  are  CLGP  FOs , In  the  case  where 
3 channels  are  shared  by  9 FOs,  segregation  of  tne 
channels  uniformly  results  in  3 FOs  sharing  a single 
channel . 

This  situation  would  eliminate  the  need  for  a 
CLGP  FO  to  change  frequency  and  would  permit  convenient, 
direct  communications  between  an  FO  and  a particular 
battery.  However  the  loss  of  flexibility  created  by 
this  arrangement  and  the  poorer  service  given  FOs 
relative  to  an  aggregated  system  having  the  same  MTTS 
does  not  appear  to  justify  the  simplification. 

The  definition  of  states  and  the  state  diagram 
for  this  case  is  snow.i  below  in  Figure  10. 


Figure  10a. 
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Figure  10b.  State  Trcnaition  Diagram 


Using  the  state  diagram  and  the  methods  outlined  on 
p.  7 one  can  write  the  state  equations  for  the  case  of 
single-channel  pre-emptive  queuing  with  3 uos . Thp^e  ai’o 
as  follows. 

State  Equations  for  the  Single-Ohannel  Case 
= - 3 ^ Pi  + a 

P2  = 3 Pf  - (u4-2\)  ^ P^ 

P3  = a P^  - (a+2\)  P^ 

P.  = 2 \ Po  - P;  -J.  P.7 

42  i\,  ( 

= 2 'v  P ^ + -i  l-'i^  - (^+l•)  P^ 

• 

r _ V I,  - * p, 

6 li.  • o 

Py  = Pg  + ^ P<,  " ^7 
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States  where  CLGP  FO  is  waiting:  l^.,  5»  6,  7 


P{CLGP  FO  is  waiting}  = + P^(t) 


+ P^(t)  + P^(t) 


Let  tje  the  expected  time  a CLGP  FO  must  wait  given 

the  system  is  in  the  k th  state  vrhen  he  arrives.  Then 


E[X^] 


■ 

k • t' 


E[X^] 


E[X^3 


i;  a-^ 


E[X^] 


Thus  the  expected  time  to  wait  given  requirement  to 
wait  is 

ix“-^  [2  P2  + P3  + U 4-  3 P^j  / ) Pk  ( 

k-2 


Parameters 

The  state  equations  for  this  case  (8)  were  solved 
for  several  values  of  u.  and  with  two  functional  foxmis 
for  ^(t).  If  one  is  concerned  about  tlie  queuing  of 
fire  requests  generated  by  outposted  Gi,GP  FOs , it  is 
necessai’y  to  considei*  a functional  form  for  \(t)  which 
is  at  least  similar  to  what  may  be  anticipated  in  a 
topical  CLGP  scenario.  It  is  anticipated  that  in  such 
a scenario,  CLGP  FOs  would  be  present  with  a company- 
sized  maneuver  element  used  as  a blocking  or  screening 
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force.  In  this  instance  the  forces  (and  the  accompany- 
ing FOs ) do  not  intend  to  hold  their  position  very  long. 

For  this  reason  the  fire  request  rate  which  will  build 
up  gradually  as  targets  begin  to  appear,  will  be  expected 
to  diminish  rapidly  at  some  point  in  time  when  the  POs 
are  compelled  to  displace  for  their  own  security.  For  tank 
targets  moving  toward  the  outpost  positions  in  assault 


T.tn  rsn  4*  ^ 


countering  mines  an  average  speed  of  advance  of  about 
5 ra/sec  is  expected.  It  is  also  expected  that  FOs  would 
withdraw  when  advancing  tanks  come  within  about  1 km 
of  their  position.  In  this  case  outposted  POs  will 
have  about  10  minutes  of  activity  before  displacing. 

The  rate  of  fire  requests  at  peak  demand  is  limited 
by  several  considerations  other  than  just  the  rate  of 
appearance  of  targets.  For  example,  the  activities 
required  of  the  PO  party,  viz.,  target  detection  and 
acquisition,  FDC  cuing,  target  cracking  and  laser  desig- 
nation also  may  limit  the  maximum  value  of  \.  Even  if 
one  member  of  the  FO  party  was  committed  to  detection 
and  acquisition  of  targets,  it  is  doubtful  whether  the 
peak  rate  for  fire  requests  v;ould  exceed  1/minute.  On 
the  basis  of  a peak  request  rate  of  l/60  (sec)"^'  achieved 
after  600  sec.,  the  following  arrival  rate  function  was 
chosen . 


\ = (1/60)  sin  Y2§0’ 


for  0 5;  t e 600  sec 


K - 0 


for  t > 600  sec 
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In  addition  to  \(t)  given  by  (12),  the  differential 
equations  (8)  were  solved  for 

\ . 1/60  , 0 s t s 600  sec 

r 0 , t > 600  sec.  (13) 

The  initial  conditions  for  all  cases  examined  are 

P;l(0)  = 1 

Pj^(O)  =0  , k 1 (lU) 


Results 

Using  \(t)  as  in  (12),  the  probability  that  a CLGP 
FO  waits  is  shown  in  Figure  11  for  several  values  of  the 
MTTS . Figures  12  and  13  were  developed  using  equation 
(13),  with  constant.  Figure  12  displays  the  probability 
that  a CLGP  FO  is  waiting  for  a channel  lo:-  the  same  set 
of  values  of  MTTS  used  in  Figure  11.  Further,  for  other 
customers  than  CLGP  FOs  who  may  wish  to  use  the  commo 
channel,  the  probability  of  not  being  blocked  is  shown 
in  Figure  13. 

These  results  clearly  indicate  that  even  for  a MTTS 
of  25  sec  (>  = 0.08  sec”^) , a single  commo  channel  is 
hopelessly  inadequate  to  satisfy  the  demands  of  three 
CLGP  FOs.  Additionally,  essentially  no  channel  capacity 
remains  for  other  than  CLGP  FOs . 
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A Commun i c aL ions  oimuiation 

A simulation  mcael  of  an  artillery  commo  system 
was  developed  to  permit  the  examination  of  a variety 
of  communications  svstem  conf  Lk’urations  with  j^reater 
ease  ana  I’lexibility  tnan  tnat  achieved  thru  analytic 
models.  The  resu,ts  of  this  simu;ation  were  compared 
with  anniv*^’r  f forv  .5 « r-r  i p rvstem  conf  1 f’l.'ra- 

tions  such  as 

a . 1 chfinr.e  1 , 3 •'  ^ 

b.  ? channels,  3 •"  2 

c.  3 channels,  ■ 

d.  3 channels,  9 r’Cs, 

each  conf i ’’urat ion  havini'  two  stares  of  service  per 
channel  and  with  only  one  type  of  message  being  trans- 
mitted. 

The  comparison  facilitated  testing  and  debugging 
the  s imu  1st '■  on  . Aside  from  the  expected  ctochaGtic 
fluetuntions  in  simulation  outputs,  there  is  essen- 
tially perfect  agreement  between  the  simulatior,  and 
analytical  results. 

The  VO  Alert 

To  appropriately  model  the  communications  traffic 
for  the  Cannon  Launched  Guided  Frojectile  System  (CLGP), 
it  is  necessary  to  treat  the  messages  transmitted  to 
the  FOs  by  the  FDC  alerting  them  to  activate  their  lasers. 
It  will  be  necessary  to  take  special  measures  to  insure 
that  the  alert  (and  its  coni'i  rmation)  got  to  the  FO  at 
the  proper  time.  However,  in  the  simulation,  it  is  not 
assumed  that  a fire  request  requires  a dedicated  channel 
thru  the  alert. 
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For  reasons  of  FO  security  it  is  desirable  to  mini- 
mize the  designation  interval.  Generally  this  designation 
interval  will  be  10  to  l5  seconds  long  as  required  by  the 
optimal  time  for  guidance.  Consequently,  a 5 second  or 
greater  lag  in  transmitting  the  alert  will  begin  to  sub- 
stantially affect  guidance  accuracy. 

To  insure  that  the  alert  is  punctual,  the  FDC  must 
allocate  a slack  interval  prior  to  tne  actual  alert  so 
as  to  insure  a timely  transmission.  This  slack  period 
blocks  incoming  message  requests  from  other  FOs  for 
the  duration  of  the  alert.  The  greater  this  slack 
period  the  greater  is  the  confidence  that  the  alert 
will  be  no  greater  than,  say,  5 seconds  late.  Because 
the  expense  of  a guidance  failure  is  great,  it  is 
anticipated  that  a high  probability  of  timely  alert 
will  be  required.  This  will  require  a non-negligible 
alert  interval  which  subtracts  thereby  from  channel 
capacity.  Clearly  a Lradeoff  must  be  made  between 
loss  of  channel  capacity  suffered  during  an  alert  and 
loss  of  guidance  accuracy  suffered  when  an  alert  is 
delayed.  Of  course,  it  is  possible  to  increase  channel 
capacity  by  formatting  messages  so  as  to  require  less 
time  (and  so  as  to  permit  interruption  by  the  FDC). 

In  the  parlance  of  queuing  theory,  this  is  equivalent 
to  the  reduction  of  service  time.  V/ith  reduced  service 
time  a higher  probability  that  the  alert  will  be  punctual 
is  achieved.  Another  element  in  this  tradeoff  is  the 
possible  addition  of  another  channel  or  dedication  of 
another  channel  to  CLGP  FOs . 

The  simulation  facilitates  tradeoffs  between  system 
parameters.  The  manner  in  which  the  alert  is  treated 
is  as  follows.  An  alert  event  i scheduled  to  be  sent 
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to  a particular  FO  on  the  channel  he  used  in  his  request 
This  alert  arrives  after  a prescribed  artillery  system 
response  time  (60  sec)  and  lasts  for  a given  alert  inter 
val.  If  the  channel  is  busy  at  this  tine,  the  alert  is 
delayed  until  that  stage  of  service  is  completed.  The 
alert  interval  was  treated  parametrically.  Results  were 
obtained  for  intervals  of  5 and  10  seconds.  These  are 
shc’vn  in  '^eble  1. 

Queuinpr  Simulation  Results 

Examples  of  queuing  simulation  results  are  shown 
in  I’igures  II4.  thru  I6.  Figure  14  shows  the  cumulative 
distribution  function  for  the  time  interval  over  which 
the  PDC  is  required  to  wait  to  transmit  an  alert. 

Figure  l5  compares  certain  statistics  from  the  analyti- 
cal model,  which  does  not  treat  an  alert,  with  a simu- 
lation using  the  alert  logic.  Figure  16  compares  the 
expected  number  of  mesearos  awaiting  transmission  for 
an  aggregated  conmo  system  of  9 FOs  and  3 channels 
with  a MTTS  of  50  sec  with  that  for  a MTTh  of  20  sec. 

All  of  the  simulation  results  were  obtained  for  the 
case  of  constant  arrival  rate  K,  starting  with  an 
empty  system. 


Table  1 

Table  of  Queuing  Sirn”.lation  Results  for  MTBA  of  60  Seconds 
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Intcjrpretation  of  the  Simulation  Results 

In  order  to  satisfactorily  handle  the  expected 
message  traffic  from  three  CLGP  FOs , it  appears  that 
two  (or  more)  channels  are  required.  A single  channel 
would  not  be  adequate  even  with  the  short  MTTS  of  20  sec. 
The  probability  that  a timely  alert  occurs  is  unaccep- 
tably small  for  a single-channel  system. 

Even  with  a two-channel  commo  system,  the  proba- 
bility of  a timely  alert  is  less  than,  say,  unless 

the  MTTS  is  20  sec  or  less  and  unless  an  adequate 
slack  period  or  alert  interval  is  allowed  prior  to  the 
latest  time  at  which  the  FO  must  turn  on  his  laser.  For 
exmiiple,  from  Table  1 we  note  that  with  a MTTS  of  20 
sec  and  a 10  sec  slack  interval,  the  probability  of  the 
FDG  being  no  more  than  5 sec  late  is  about  91^  even  tho 
the  expected  number  in  the  queue  is  only  0.05« 
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a.  MER,  DRSAR-SAM,  23  Jan  76,  subject:  Description  of  Computer 

Program:  Exterior  Ballistics  of  Boosted  Rockets.  This  MER  appears 

in  the  Systems  Analysis  Directorate  Activities  Summary  January  - 
February  1976,  DRSAR/SA/N-42. 

b.  Letter,  SARRI-LR  to  AMXBR-EB-FT,  8 Jan  76,  subject:  CLGP 

Terminal  Homing  Models. 

c.  MER,  DRSAR-SAM,  20  Apr  77,  subject:  An  Imprc”cd  Algorithm  for 

the  Glide  Mode  of  Copperhead  Used  in  a 3 DOF  Flight  Simulation  (Incl  1). 

2.  The  modified  three-degree-of-freedom  (3  DOF)  ballistic  simulation 
described  in  Ref  la  was  sent  to  BRL  in  January  1976.  The  computer 
program  was  intended  to  be  an  efficient  means  for  generating  range 
tables  for  Copperhead,  including  the  glide  (FUFO)  mode.  In  achieving 
this  objective  a very  simple  al gori thm--gl ide-hold , discussed  in  Ref  la-- 
was  used  to  simulate  the  FUFO  mode  of  flight. 

3.  As  delineated  in  Ref  lb,  good  agreement  was  obtained  between  ranges 
generated  by  the  ARMCOM  Systems  Analysis  Model  (SA  3)  and  the  elaborate 
Rodman  six-degree-of-freedom  model  (ROD  6).  However,  it  is  believed 
that  the  agreement  for  glide  trajectories  between  SA  3 and  ROD  6 is 
somewhat  misleading  since  all  glide-mode  comparisons  were  made  for 
firiigs  in  the  lower  register  and  for  shots  in  which  the  glide  angle 
did  not  change  greatly  from  the  nominal,  twenty-degree  value. 

4.  In  order  to  improve  agreement  between  these  models  for  off-nominal 
conditions  using  the  glide  mode,  SA  3 has  been  changed  to  mechanize  a 
new  glide  algorithm.  This  algorithm  is  described  in  Ref  Ic  (Incl  1) 


49 


DRSAR-SA 

SUBJECT:  Computer  Simulations  for  the  Ballistics  of  Copperhead 


where  results  are  compared  with  results  of  6 DOF  models  used  by  Martin 
Marietta  Corporation  (MMC  6)  and  by  Rodman  Labs  (ROD  6). 


5.  Our  computational  experience  with  this  version  of  SA  3 (April  1977) 
indicates  that  trajectory  shape,  particularly  glide  angle,  better  represents 
results  of  6 DOF  models  than  the  1976  version.  To  expedite  your  receipt 
of  this  model,  the  source  program  card  deck  with  listing  is  being  sent 
under  separate  cover. 

FOR  THE  COMMANDER: 

EJGNCD 

1 Incl 
as 

OTTO  F.  HAASE,  JR. 

Acting  Director 

Systems' Analysis  Directorate 

CF: 

Project  Manager,  Cannon  Artillery  Weapons  Systems,  ATTN:  DRCPM-CAWS/ 

Mr.  E.  Manley/Mr.  E.  Zimpo/COL  R.  Nulk,  Dover,  NJ  07801 
Project  Manager,  Cannon  Artillery  Weapons  Systems,  ATTN:  DRCPM-CAWS-SIA/ 

Mr.  J,  Williams,  Dover,  NJ  07801 

Commander,  White  Sands  Missile  Range,  ATTN:  STEWS-TE-PC/Mr.  S.  Kadner, 

Bldg  1544,  White  Sanas,  NM  88002 

Commander,  US  Army  Armament  Materiel  Readiness  Command,  ATTN:  DRDAR-LCA-PP/ 

Dr.  Amoruso,  Rock  Island,  IL  61201 
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MEMORANDUM  FOR  RECORD 

SUBJECT:  An  Improved  Algorithm  for  the  Glide  Mode  of  Copperhead  Used  in 

a 3 DOF  Flight  Simulation 


1.  References: 


a.  MFR,  DRSAR-S7U-I,  23  Jan  76,  subject:  Description  of  Computer 

Program:  Exterior  Ballistics  of  Boosted  Rockets.  This  MFR  appears  in 

the  Systems  Analysis  Directorate  Activities  Summary  January  - February 
1976,  DRSAR/SA/N-42. 


b.  Letter,  SARRI-LR  to  AMXBR-EB-FT,  8 Jan  76,  subject: 
Homing  Models. 


CLGP  TerLiinai 


2.  Previous  Work. 

In  modeiing  the  glide,  iUFO,  or  attitude-hold  mode  of  Copperhead  witl.in 
the  general  purpose,  point-mass,  three-degree-cf-f reedom  (3  DOF)  program 
EXBAL,  one  generaliy  desires  computational  efficiency  "’nd  fidelity  of 
pertinent  flight  variables  to  the  values  obtained  from  more  elaborate 
six-degree-of-freedom  (6  DOF)  models.  In  Ref  la  an  algorithm  is  presented 
in  which  the  glide  angle  is  fixed  or  held  by  appropriate  choice  of  aero- 
dynamic forces.  Two  options  are  permitted.  In  one,  the  glide  angle  is 
held  after  a prescribed  glide  angle  is  reached.  In  the  other  option  the 
glide  angle  achieved  after  a prescribed  time,  TENABL,  is  held  for  the 
remainder  of  the  flight. 


3.  Aitho  tlie  glide  mode  of  Copperhead  does  not  actually  operate  in  this 
manner.  Ref  lb  shows  that  range  estimates  of  good  accuracy  can  be  obtained 
using  Lliis  algorithm.  However,  several  shortcomings  of  this  algoritkmi 
have  been  observed.  For  one,  the  glide  angle  actually  assumed  in  the  glide 
mode  In  not  constant  but  takes  cn  a value  such  as  to  preserve  body  attitude. 
Parenthetically,  note  that  at titu Je-liold  is  most  descrjptive  of  this  mode 

of  flight.  Thus,  the  glide  angle  is  only  appicximatcly  constant  and  may 
depart  .significantly  from  the  nominal  -20  degrees  de.sircd.  For  these 
trajectories,  a glide-hold  aigorit’nm  does  not  adeqtiately  represent  trajectory 
chape. 

4,  .brother  shortcoming  of  the  earlier  ( 976)  glide  algorithm  is  the  absence 
of  n direct  input  of  the  Copperiu.ad  timi..  setting  used  by  the  6 DOF  models. 
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Tq.  This  timer  setting  represents  the  time  from  launch  at  which  the 
30  V main  projectile  power  becomes  available.  The  time  Tq  is  the  (integer) 
time  (in  seconds)  set  by  code  switch  into  the  projectile  plus  the  approximate 
0.5  second  interval  to  account  for  the  following  lags  not  included  in  the 
projectile  timer:  (a)  the  approximately  0.30  second  interval  for  the  11  v 

battery  to  activate  loliowing  launch  and  the  approximately  0.20  second 
interval  required  for  the  30  v battery  to  achieve  acceptable  power  after 
it  is  activated.  All  ether  events  occurring  within  the  projectile  are 
sequenced  from  the  base  time  Tq.  These  events  are  tabulated  below: 

SEQUENCE  OF  EVENTS  IN  COPPERHEAD  TIMER 


Event 

No. 

Time 

(s) 

Description 

1. 

To 

timer  sequence  starts  (main  power  at  + 30  v) 

2. 

Tq+2 

roll  control  starts 

3. 

Tq+3 

initiate  seeker  gyro  spinup 

4. 

Tq+4 

attitude-hold  enable  switch  is  closed 

5. 

Tq+5 

start  to  extend  wings,  apply  g bias,  and  free 
gyro 

5.  Altho  internal  degrees  of  freedom  within  the  projectile  are  not  simulated 
with  the  3 DOF  model,  it  is  important  to  properly  incorporate  their  kineiLStic 
effect.  In  contrast  to  the  earlier  algorithm,  one  option  of  the  improved 
algorithm  described  below  treats  the  events  which  occur  after  Tq  in  a manner 
similar  to  that  of  the  actual  projectile.  As  with  the  previous  algorithm, 
the  principal  assumption  employed  is  that  the  controlled  projectile  always 
remains  in  a trim  condition. 

6 . Options  for  the  Imp rove d Program . 

There  are  three  sv/itcViable  options  associated  with  the  glide  mode 
in  eurrent  version  of  EXB^L.  All  opti'^ns  require  the  glide  mede  switch 
IGLIDE  to  be  set  to  unity.  In  option  1 a desired  glide  angle  (GLIDE) 
is  provided  and  the  program  calculates  t e body  angle  of  attack  required  to 
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hold  body  attitude  after  reaching  the  desired  glide  angle.  Before  attitude 
hold  of  the  body  is  initiated  an  initial  angle  of  attack  is  calculated 
such  that  the  lift  force  equals  the  gravitational  force  normal  to  the 
velocity  vector  at  that  point  in  the  flight.  See  Attachment  1 for  details. 

This  option  may  be  useful  when  timer  settings  are  unavailable. 

7.  In  option  2 a time  TENABL  (in  secs)  is  provided  and  CLIDE  is  set 
to  -90  degrees.  With  this  option  the  calculations  for  angle  of  attack 
to  hold  body  attitude  start  at  TEFABL.  A value  of  TEUABu  equal  to  tiiat 
at  which  the  velocity  angle  0 = GLIDE  produces  a trajectory  identical 
to  that  of  option  1. 

8.  Option  3 requires  the  program  user  to  provide  time  Tq,  called  TMO'LMC 
in  t'ne  program.  This  option  is  exercised  by  setting  the  sv.’itch  15MCSW 

to  unity.  So  far  as  projectile  kinematics  are  concerned,  the  projectile 
remains  ballistic  until  event  number  4.  When  event  number  4 occurs  at 
To+4(s),  a control  surface  deflection  in  pitch,  6p,  is  calculated  con- 
tinuously and  applied  thruout  the  interval  betwecii  events  4 and  5.  This 
algorithm  bases  the  calculation  of  6p  upon  the  equilibrium  behavior  of 
the  projectile  and  autopilot  between  events  4 and  5.  This  behavior  is 
approxihiatsly  one  in  which  th.c  pitch  .ic.fJcctior.  of  the  fins  end  the  nscoc  iat-f'ri 
trim  angle  of  attack  are  proportional  to  the  average  turning  race  of  the 
• 

velocity  vector,  0.  Thus, 


6 = K 0 (K,  a constant) . 

P 

K is  approx.  0.5  deg/ (deg/sec) . The  turning  rate  0 is  calculated  continuously 
by 

• •••  •••  2 

0 = (x  y - y x)/v  . 

An  average  value  of  0 is  calculated  by  exponential  smootliing.  This,  of 
course,  assumes  that  the  autopilot  and  projectile  system  er-dilbit  appro;-.imately 
first-order  dynamics  between  the  events  4 and  5.  Specifically, 


with  T a 
value  of 
equation 


time  constant  and  s the  Laplace  differential  operator.  The 
T is  apprc'ximatel y 0.2  secs.  Using  exponential  smcociiing,  this 
is  solved  implicitly  by  the  recursive  procedure: 
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0 = c 6,  , + C-0. 

i o i-1  1 1 


with 


-h/T 


c = e 
o 


Cl  = 1 - Co  . 


where  h is  the  integration  time  step,  ie, 


t.  = t.  , + h 
1 1-1 


Having  calculated  6p,  the  trim  angle  of  attack  («(■)  is  obtained  by  linear 
interpolation  in  the  following  Cable. 


TRIM  ATTACK*  VERSUS  MA.CII  NUMBER  AND  CONTROL  DEFLECTION 
Entries  are  Oj.  (deg) 


Mach 

Number 

6 

(deg) 

0 

5 

10 

15 

0.5 

0 

7.4 

12.1 

17.7 

0.8 

0 

7.2 

11.8 

17.2 

0.9 

0 

7.4 

11.6 

17.1 

1.0 

0 

1 

7.1 

1 

11.3 

17.1 

* Reference:  Wind  Tun:icl  Data  Analysis  of  3/4  Scale  Model  for  the 

XM712  Projectile,  27  Feb  76,  pp.  72-75. 
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9.  With  all  program  optiohs  the  normal  body  force  is  calculated  from 
the  normal  force  coefficient  at  trim.  This  coefficient  is  a function 
of  Mach  number,  M,  and  trim  attack,  a^,  ie. 


Values  of  C^  are  obtained  by  two-way  linear  interpolation  in  a table  as 
explained  in  Attachment  1. 

10.  Comparison  of  EXBAL  Results  With  6 DOF  Models. 

Several  different  Copperhead  trajectories  have  been  generated  using 
the  improved  glide-mode  algorithm  in  SA  3 for  the  purpose  of  comparirig 
these  results  with  those  of  the  6 DOF  models.  Results  of  three  glide 
trajectories  and  a ballistic  trajectory  are  shown  in  Table  3.  Comparable 
results  for  the  Rodman  (ROD  6)  and  Martin  Marietta  Corporation  (MMC) 

6 DOF  models  are  shown,  where  available.  For  all  glide  trajectories 
with  SA  3,  option  number  3 was  used.  The  ballistic  option  is  compared 
with  results  of  ROD  6 to  demonstrate  the  virtual  equivalence  of  program 
results  when  using  this  option. 

11.  Altho  the  range  generated  by  SA  3 is  quite  comparable  with  that  of 
ROD  6 and  MMC  ^ glide , it  is  noted  that  the  final  flight  speed  is  about 
1 to  2 m./s  less  than  the  6 DOF  results  and  the  final  glide  velocity  angle 
is  about  0.1  to  0.3  degree  less.  For  the  66.5  deg  QE  shot  a trajectory 
is  displayed  in  Figure  1.  The  range  given  by  SA  3 is  somewhat  longer 
and  the  apogee  higher  than  that  of  ROD  6.  This  in  part,  accounts 

for  the  0.4  sec  longer  time  of  flight  for  SA  3. 

12.  A significant  flight  variable  for  the  glide  mode  is  the  glide  angle,  0. 
This  is  shown  versus  time  for  both  SA  3 and  ROD  6 in  Figure  2.  During 

the  transition  interval  from  ballistic  flight  to  attitude-hold  flight — 

44.5  (s)  to  about  46  (s) — there  is  excellent  agreement  between  these  models. 
At  the  end  of  the  flight  the  SA  3 result  becomes  about  0.3  degree  less 
(steeper) . 

13.  Conclusions . 

In  view  of  the  excellent  agreement  displayed  here,  one  is  led  to 
believe  that  the  improved  algorithm  for  the  glide  mode  will  prove  suit  iOle 
for  range  table  generation.  Further,  the  present  mechanization  does  not 
require  substantially  more  CPU  time  than  the  earlier  (1976)  version. 
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In  fact  this  program  is  still  one-to-two  orders  of  magnitude  faster  than 
the  6 DOF  models.  . This  suggests  that  it  could  become  the  basis  for  a real- 
time F/C  algorithm. 


1 Incl 

• Attachment  1 


< i-  jZ  L .•  t 

GEORGE  J.  SCHLENKER 
Operations  Research  Analyst 
Methodology  Division 
Systems  Analysis  Directorate 
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COMPARISONS  OF  Rf:SULTS  FROM  COPPERHEAD  SIMULATIONS 
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yig-irc  1.  Reference  Trajectory  for  Copperl'ead  in  Glide  Mode  Using  the  SA  3 Flight  Simulation 


ATTACKiMlLNT  1 


Al^^orlthm  for  Att itude-Hold  Logic  In  Copp e rhead  3D0F  Simulation 

For  simulations  in  which  a preset-time  option  is  not  used,  the 
program  determines  the  time  for  commencement  of  attitude  hold.  Initially 
the  glide  angle  0 will  approach  the  desired  glide  angle  6^  a]gebralcally 
from  above.  (0^  is  negative)  l\Tien  9 - Gq  £ ^ positive),  an  initial 
angle  of  attack,  a^,  is  computed  such  that  the  associated  lift  component 
equals  the  component  of  gravity  normal  to  the  velocity  vector,  ie,  such 
that 

F-  = li  g cos  6 
L p ^ 

with 


and 


Ft  = F„  cos  a 
L N o 


F„  - 0,(0)  A„ 


with  projectile  mass  and  reference  area  A and  dynamic  pressure  q. 

Since  is  small — typically  less  than  10° — an  iterative  jirocedure 

is  employed  in  which  the  first  iteration  assumes  that  = F„  so  that 

L N 

(1)  M 8 cos  0 

C„(0)  = 

N Aq 

The  value  of  is  obtained  by  Interpolation  in  the  tabular  function 

Cj^.(H,  a^)  with  the  trim  angle  of  attack. 


Thus , 


* (1) 

C^(0)  = C^.(M*, 


where  M*  is  the  local  Mach  number.  The  interpolation  procedure  is  described 
below. 

Then,  form  the  second  iterate  for  lIiu  normal  fcrce  coefficient: 

(2)  M 8 cos  0 

" AO 

Aq  cos  a 
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(2) 

The  value  of  obtained  by  requiring  that  = ^^^(0)  is  taken 
as  the  initial  trim  angle  of  attack  at  the  start  of  attitude  hold. 

Interpolation  Procedure 

To  calculate.  , interpolate  on  Mach  number  in  the  a^)  table 

obtaining  C^(M*,  a^)  at  the  local  Mach  number  M*  for  values  of 

= {O,  5,  10,  15  (deg)}.  Then,  is  obtained  by  linearly  interpolating 
with  0^(0)  as  argument  in  this  table. 

The  value  of  initial  body  attitude  is  given  by 

\<») 

For  subsequent  calculations  during  the  attitude-hold  trajectory, 

the  value  of  angle  of  attack,  a,  is  calculated  which  preserves  body 

attitude,  ie , for  v;hich  0.  = 0.(0). 

b b 

Thus , 

a = 0,.(0)  - 0. 

D 

If  the  above  value  of  a satisfies 

a < a^^(M) , the  maximum  trim  angle  at  M '■he  instantaneous  value 
of  Mach  number,  the  lift  force  is  computed  as 

F = C (a,  M)  Aq  cos  a 

with  obtained  by  two-way  Interpolation.  Otherwise,  a is  limited  to 

(M)  and  lift  is  calculated  as 
tm 

F = C,,  (a  , M)  Aq  cos  a 
L N'  tm’  ^ tni 

Induced  drag  due  to  lift  is  calculated  as 

Fpj.  = - C^,(a,  M)  Aq  sin  a . 
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TABLE  1.  MAXIMUM  TRIM  ANGLE  OF  ATTACK 
VERSUS  MACH  NUMBER  (AT  6 = 12°) 


Mach 

Number 

(max) 

(deg) 

0.5 

14.3 

0.8 

14.0 

0.9 

13.8 

1.0 

13.6 

TABLE  2.  NORMAL  FORCE  COEFFICIENT 
VERSUS  MACH  NUMBER  AT  TRT’.;  ATTACK 


Cj/M,  a ) 


DRSAR-CG 


DRSAR-SA 


2 ; 1£// 

H.W.  Michels/jls/5891 


Quality  Levels  and  Effectiveness 


To  inform  the  CG,  ARRCOM,  of  examples  of  commodity  requirements  which,  in  analysis, 
provide  only  marginal  benefit  in  operational  cost-effectiveness  at  the  force  level. 


1.  Reference  Project  List  Review,  MG  Eicher  and  Mr.  Rhian,  DRSAR-SA,  28  March  1977; 
question  on  Quality  Control. 

2.  We  have  evaluated  the  influences  of  system  alternatives  on  effectiveness  achieved, 
mix  of  ammunition  expended,  and  its  dollar  cost.  We  used  our  BLUE  artillery  force  to 
engage  total  target  arrays  subject  to  TRADOC-approved  mission  scenarios,  tactics,  and 
usage  doctrine.  Results  showed  that  some  specified  commodity  requirements  (performance 
attributes/improvements)  were  not  justified  by  the  realization  of  a significant  im- 
provement in  operational  cost-effectiveness  of  the  field  artillery.  Examples  (paras 

3.  4 and  5)  are  summarized  below.  ^ 

3.  The  155mm,  HE,  M549  RAP  production  rounds  exhibited  a greater  range  probable  error 
(0.4%  of  range)  than  the  developmental  rounds  (0.3%  of  range).  Analysis  showed  that 
more  precise  (developmental)  range  variation  significantly  improved  effectiveness  only 
against  small  area  targets  engaged  by  observer-adjusted  fire.  On  the  basis  of  the 
Legal  Mix  target  arrays,  these  engagement  conditions  were  predicted  to  occur  in  only 
15  percent  of  the  M549  RAP  firings.  In  the  remaining  85  percent  of  the  firings,  when 
engaging  larger  targets  and/or  using  unadjusted  fire,  n£  significant  improvement  was 
expected  when  using  the  developmental  variation.  This  analysis  was  reported  in  a 
briefing  to  the  CG,  ARMCOM,  in  April  1974. 

4.  The  improved  precision  specified  for  two  time  fuzes  (the  XM724  ET  and  the  M577  MT), 
which  are  candidate  fuzes  for  the  second  generation  dual-purpose  ICM's,  did  not  result 
in  a significant  improvement  in  artillery  force  level  operational  effectiveness  nor  in 
a significant  ammunition  cost  savings  when  compared  with  an  alternative  in  which  the 
ICM  dispersion  was  based  on  the  use  of  the  less  precise  M565  MT  fuze.  This  analysis 
was  published  in  report  DRSAR/SA/R-20,  July  1976. 

5.  The  155mm,  HE.XM708  hi-frag  projectile  provides  signif icantly  improved  single 
round  effectiveness  on  personnel  targets  in  comparison  to  the  155mm,  HE,  M107  pro- 
jectile. However,  analysis  showea  that  the  XM708  yielded  n£  significant  effectiveness 
or  cost  advantage  when  included  in  the  ammunition  mix  at  the  artillery  force  level 
(i.e.,  in  comparison  to  the  results  with  the  M107  in  the  mix).  This  analysis  was 
published  in  technical  note  DRSAR/SA/N-49 , December  1976. 

6.  The  155mm,  HE,M483A1  ICM  (dual-purpose,  second  generation  improved  conventional 
munition)  is  an  example  of  a commodity  improvement  which,  in  analysis,  doe£  signifi- 
cantly benefit  operational  cost-effectiveness  at  the  force  level.  An  analysis  showed 
that  the  M483A1  provides  essential  anti-armor  capability  among  the  conventional, 
unguided  155mm  field  artillery  munitions.  Another  analysis  showed  that  only  one- 
eighth  as  many  battery  volleys  (fewer  for  some  conditions)  of  M483Al's  are  required  to 
achieve  the  same  level  of  effectiveness  against  area  personnel  targets  as  that  achieved 
by  battery  fire  of  155mm,  HE,  M107  projectiles.  The  analysis  pertinent  to  anti-armor 
capability  was  documented  in  a DRSAR-SAA  Memorandum,  November  1976.  The  other  analysis 
was  documented  as  a AMSAR-SAA  Memorandum,  dated  April  1975. 
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7.  The  failure  of  one  or  more  specific  requirements  (attributes  or  quality  levels) 
of  a commodity  to  be  justified  in  terms  of  improved  operational  cost-effectiveness 
does  not  necessarily  imply  that  the  commodity  should  not  be  acquired.  For  instance, 
the  justification  for  continued  development  of  the  XM724  electronic  time  fuze  as  a 
candidate  fuze  for  the  new  ICM's  (re:  paragraph  4 above)  was  based  on  safety  and 

the  general  evaluation  of  the  electronic  timing  industry 


M.  RHIAN 

Director,  Systems  Analysis  Directorate 


CF: 

DRSAR-QA,  J.J.  Obren 
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